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Abstract 

We give necessary conditions for the efficient simulation of both bipar- 
tite and multipartite Hamiltonians, which are independent of the eigen- 
values of Hamiltonians and based on the algebraic-geometric invariants 
introduced in [1] and [2]. The results show that the problem of efficient 
simulation of Hamiltonians on arbitrary bipartite or multipartite quantum 
systems cannot be described by only using eigenvalues, which is quite dif- 
ferent to the two-qubit case. 

Historically the idea of simulating Hamiltonian time evolutions was the first 
motivation for quantum computation [3] . Recently the ability of nonlocal Hamil- 
tonians to simulate one another is a popular topic , which has applications in 
quantum control theory [4] , quantum computation [5] , [6] , [7] , [8] and the task of 
generating entanglement [9] [10]. The problem to parameterize the nonlocal 
properties of interaction Hamiltonians, so as to characterize the efficiency with 
which they can be used to simulate one another, is theoretically and experimen- 
tally important. There have been very active research on this problem ([11], 
[12], [13], [14], [15], [16], [17]). For the general treatments of this topic , we refer to 
[11] as the main reference. 



In [11] it was shown that the efficiency with which Hamiltonian H , together 
with local operations , simulates another Hamiltonian H' can be used as a crite- 
rion to endow the set of Hamiltonians with a partial order structure, that allows 
to compare the nonlocal capabilities of H and H' . For two-qubit Hamiltonians, 
it was shown that the problem of simulation of Hamiltonians can be reduced 
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to the case of so-called normal forms of these Hamiltonians (Theorem in sec- 
tion V A of [11]). For these Hamiltonians in their normal forms , a beautiful 
necessary and sufficient condition about the simulating capability in terms of 
the majorization of eigenvalues of these Hamiltonians was given. This indicated 
that the partial order structure endowed on the two-qubit Hamiltonians is in 
close analogy to the partial ordering of bipartite pure states endowed by their 
capabilities to be converted by LOCC ([18]). 

It is natural to consider the simulation problem of Hamiltonians on arbi- 
trary bipartite quantum systems. We can imagine that in higher dimensions, 
Hamiltonians have more nonlocal degrees of freedom than the two qubit case, 
and there is no result about "normal forms" of Hamiltonians on arbitrary bipar- 
tite quantum systems, this may make the problem more difficult and it seems 
hopeless to give a characterization based only on the majorization of some nu- 
merical quantities (as in Theorem of section F of [11]). In [19] it is proved 
that TO^n^ — m-^ — n"^ + 1 nonlocal parameters are needed to describe the set of 
equivalent classes of bipartite mixed states on (g) Hg under local unitary op- 
erations. A similar parameter counting argument as [19] shows that there must 
be at least m^n^ — 2(m^ -|- n^) + 3 nonlocal parameters for the equivalent classes 
of bipartite Hamiltonians on H^^H^ under local unitary operations. When we 
consider the ability of nonlocal Hamiltonians to simulate one another with the 
help of local unitary operations, it is natural to imagine that these continuous 
invariants of Hamiltonians under local unitary operations may give constraints 
on these Hamiltonians H and H' if there is a simulation relation using local 
unitary operations between them ,and any such constraint must be expressed 
by these invariants. In our previous works [1] and [2] the algebraic-geometric 
invariants of bipartite mixed states (i.e. semi-positive sclf-adjoint operators) 
were introduced as their nonlocal invariants (i.e., these algebraic sets are kept 
invariant under local unitary operations). These algebraic-geometric invariants 
depends only on eigenvectors and are independent of eigenvalues of the semi- 
positive adjoint operators. We can think these algebraic-geometric invariants 
as nonlocal invariants of semi-positive bipartite Hamiltonians and ask if there 
exists any constraint on these invariants of two semi-positive Hamiltonians H 
and H' if H can be simulated by H' using local unitary operations. 

In this paper, we show that the efficient simulation relation between two 
semi-positive bipartite Hamiltonians of the same rank implies the equalities of 
these algebraic-geometric invariants of them. This necessary condition is also 
extended to the efficient simulation of multipartite Hamiltonians. Since these 
algebraic-geometric invariants are independent of the eigenvalues and only mea- 
sure the position of eigenvectors of the Hamiltonians. Thus our results strongly 
suggest that the eigenvectors play a more fundamental role in the efficient sim- 
ulation problem of Hamiltoians on arbitrary bipartite or multipartite systems. 
This is quite different to the two qubit case studied in [11]. 
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In this paper, we say, for two bipartite Hamiltonians H and H' on H™^H^, 
H' can be efficiently simulated by H with local unitary operations, write as 
H' ~<LU H, if H' can be written as a convex combination of conjugates of H by 
local unitary operations, H' = pi{Ui (g) Vi)H{Ui (g) Fi)^ + ... +Ps{Us (g) Vs)H{Us (g 
Vs)^ , where are positive real numbers such that pi + ... + pg = 1 , 

Ui,...,Us and Vi, Vg are unitary operations on iJ™ and Hg respectively. Here 
we use t for the adjoint. This is equivalent to the notion "infinitesimal simu- 
lation" in [11] and "first order simulation in [16]. In [11] and [16] it is shown 
that "local terms" like I (g) Kb and Ka (g / are irrelevant to the simulation 
problem upto the second order, thus they consider the simulation problem for 
Hamiltonians without local terms' effect. Our definition here is more restricted 
without neglecting the local terms. 

If H' -< Lij H, where is a semi-positive self-adjoint operator, it is clear that 
H' has to be a semi-positive self-adjoint operator and tr{H') = tr{H). Thus it 
is clear that H' -<lu H is equivalent to H' - ^^^^^Imn ~<lu H - ^^^Imn, so 
we do not restrict to the traceless Hamiltonians in this paper. 

We have the following observation. First we recall the following result in [20] . 

Lemma 1. Let T = T,iPi\vi){vi\, where pi's are positive real numbers, be 
a positive self-adjoint operator on a finite dimensional Hilbert space. Then the 
range range{T) of T is the linear span of vectors Vi 's. 

li H = \v){v\ and H' = \v'){v'\ where \v) and \v') are pure states and H' can 

be simulated by H efficiently, ie., H' H, actually the Schmidt ranks of \v) 
and \v') have to be the same. In fact, if there exist positive numbers pi, and 
local unitary operations ?7i(g Vi, Us<SiVs, such that, T,iPiUi^ViH{Ui0Vi)^ = 
H', it is clear that Ui'S)ViH{U.i_®Vt)'' = \{U,,^'Vi)v) {{Ut®Vi)v\, and from Lemma 
1, \{Ui (g Vi)v) is in the range of H' . Hence \v') = \ {Ui (g Vi)v) and the Schmidt 
ranks of \v) and \v') have to be the same. 

For semi-positive bipartite Hamiltonians (equivalently ,bipartite mixed states, 
i.e., semi-positive self-adjoint operators) on i?™ ® algebraic sets V^iH) in 
(jpm-i (^respectively Vg{H) in CP'^~^) are introduced in [1] as the degenerating 
locus of the measurement of them by separable pure states. For any given semi- 
positive self-adjoint operator (bipartite mixed states or semi-positive Hamilto- 
nians) p on _ff™ (g Hg , we consider the expression (g ^2|p|'^i ^ 4'2) for any 
pure states 0i e ff^ and 02 G H^- For any fixed 0i G where PiH"^) 

is the projective space of all pure states in H'^, (0i (g 02|p|0i <g ^2) is a Hermi- 
tian bilinear form on Hg, denoted by {(f>i\p\(j)i) . We consider the degenerating 
locus of this bilinear form, ie., V^(p) = {4>i S P{H^) : ranA:(((/)i |p|(/>i)) < k} 
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for = 0, 1, ...,n — 1. We can use the coordinate form of this formahsm. Let 
{|11), |ln), |ml), |mn)} be the standard orthogonal basis of (g) Hg 
and p be an arbitrary (scmi)positive self-adjoint operator. Wc represent the 
matrix of p in the basis {|11), ...|ln), |ml), Imn)}, and consider p as a 
blocked matrix p = {pij)i<i<m,i<j<m with each block pij a. n x n matrix cor- 
responding to the \in) rows and the \jn) columns. For any pure 
state (pi = ri|l) + ... + rm\ni) G P{H^) the matrix of the Hermitian hnear form 
(^i|p|0i) with the basis |1), \n) is TiijnrjPij. Thus the "degenerating locus" 
is actually as follows. 

V^ip) = {(n, ...,rm) e (7P™-i : rank{J:ijrir]pij) < k} 

for A; = 0, 1, ...,n — 1. Similarly Vg{p) C CP"~^ can be defined. It is known 
from Theorem 1 and 2 of [1] that these sets are algebraic sets (zero locus of 
several multi- variable polynomials, see [21]) and they are invariants under local 
unitary operations depending only on the eigenvectors of p. Actually these al- 
gebraic sets can be computed easily as follows. 

Let {|11) , .... , .... |ml), |mn)} be the standard orthogonal basis of 
as above and p = ^l-iPi\vi) {vi \ be any given representation of p as a convex 
combination of projections withpi, ...,pt > (for example, we can take the spec- 
tral decomposition p = T,l^iXi\ipi) {ipi\ as a such representation). Suppose vi = 
^^jZiCiijilij) , X = {aiji)i<i<m,i<j<n,i<l<t is the mnxt matrix. Then it is clear 
that the matrix representation of p with the basis {|11), |ln), |ml), \mn)} 
is XPX^ , where P is the diagonal matrix with diagonal entries pi, ...TPt- We 
may consider the mn x t matrix X as a m x 1 blocked matrix with each block X^, 
where m; = 1, m, a nx t matrix corresponding to \wn)}. It is clear 

Pij — XiPXj and T,ijrir^jPij = {'EiriXi)P{T,iriXi)^ . From simple linear alge- 
bra V^{p) is just the set of points (ri, r^) in CP"^~^ such that rank{YjiriXi) 
is less than fc -|- 1, ie., V^{p) is the algebraic set in CP™"^ as the zero locus of 
the determinants of all (A; -|- 1) x (fc -|- 1) submatrices of SjrjXi (see [1]). 

We can sec that V^{p) is independent of eigenvalues Ai's if it is computed 
from the spectral decomposition p = S[^]^Ai|'!/'i)(V'i|; since it is computed from 
the matrix X depending only on eigenvectors |V^i)'s. For example, let T be a 
two-qubit mixed state with some of the following 4 Bell states as its eigenvectors. 

ki) = 7l(|ll) + |22)) 
h) = ^(|ll)-|22)) 
|t;3) = ^(|12)H-|21)) 
K) = 45(|12)-|21)) 
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It is easy to calculate the matrix nXi + r2X2 of the whole 4 Bell states, it 
is the following 2x4 matrix. 



Therefore V^{T) and V^{T) can be computed from the submatrix consisting 
of rank{T) columns of the above matrix. Thus V^{T) is always empty, V^{T) 
is the set of 2 points when rank{T) = 2 and empty when rank{T) = 3, 4. 

Prom [1], Schmidt ranks of pure states p(ie, projection operator to a unit vec- 
tor) are just the codimcnsions of the algebraic sets {codimV^{p) = codimVg^fj)) . 
Therefore it is natural to think the above observation can be extended to the 
equalities of these algebraic sets of arbitrary bipartite semi-positive Hamilto- 
nians of the same rank if they can be simulated efficiently. In this paper we 
give such a necessary condition about the efScient simulation of semi-positive 
Hamiltonians. 

Theorem 1. Let H and H' be the semi-positive Hamiltonians on the bi- 
partite quantum system (g) with the same rank, ie., dim{range{H)) = 
dim{range{H')). Suppose that H' ^lu H, that is , H' can be simulated by 
H efficiently by using local unitary operations. Then V^{H) = V^{H') for 
k = 0, n — 1 and Vg{H) = Vg(H') for k = 0, m — 1, here the equality of 
algebraic sets means they are isomorphic via projective linear transformations 
of complex projective spaces. 

The following observation is the the key point of the proof of Theorem 1. 
From Lemma 1 in [20] as cited above, the range of p is the linear span of vectors 
\vi), \vt). We take any dim{range{p)) linear independent vectors in the set 
{\vi), \vt)}, say they are \vi), \vs) , where s = dim{range{p)) . Let X' be 
the mn x s matrix with columns corresponding to the s vectors \vi), ...,\Vs)'s 
coordinates in the standard basis of iJ™ (8) Hg. Then X' is a submatrix of 
the above-described matrix X and each column of X is a linear combination of 
columns in X' . We consider X' as m x 1 blocked matrix with blocks X[, X'^ 
{n X s matrices) as above. It is clear that V^{p) is just the zero locus of de- 
terminants of all (fc + 1) X (k + 1) submatriccs of E^riX-, since any column in 
T,iriXi is a linear combination of columns in SjrjX- ( thus rank(Y^iriXi) < k 
is equivalent to rank{T,iriX[) < k). 

Proof of Theorem 1. Suppose H' -<lu H, then there exist positive 
numbers pi, ■■■,Ps and local unitary operations Ui (g) Vi, Ut (8) Vt, such that, 

^^iPrU.,®V,H{Ui»V^y = H'. Let H = T.l^^q,\tlj,){iP,\, where s = dim,{range{H)) 
, qi,...,qs are eigenvalues of H and | V'l ),••-, IV's) a-re eigenvectors of H. Then 
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it is clear that ([/» (g) V^)H{U^ ® Vi)^ = ^^^qj\{U^ ® Vi)il)j){{U^ ® Vi)il)j\ and 
thus H' = T,l"^i .j^iPiqj\{Ui (8) Vi)ij)j){{Ui ® Vi)il)j\. This is a representation of 
H' as a convex combination of projections. From Lemma 1 range{H') is the 
hnear span of \{U\®V\)'il)\) , |(C/i (g) Vi)'!/'s) since they are dim{range{H'j) = s 
hnear independent vectors in range{H'). Prom our above observation V^(iJ') 
can be computed from the matrix X' of vectors \{Ui (E) Vi)tpi), \{Ui Vi)tps) 
and thus V^{H') = V^{{Ui ® Vi)H{Ui O Vi^) from the definition. Thus the 
conclusion follows from Theorem 1 in [1] . 

Since the algebraic-geometric invariants are independent of eigenvalues, thus 
our above theorem is a necessary condition of simulation of Hamiltonians with- 
out referring to eigenvalues. As described in e.g. [11] and [16], local terms like 
I ® Kb and Ka ® I are considered irrelevant to the simulation process upto 
the second order, this leads to the so-called normal forms of two-qubit Hamil- 
tonians. Wc can recall the Theorem in section F of [11], for Hamiltonians H 
and H' in their normal forms, ie., H = HihiOi ® Oi and H' = Ti^h'^Gi ® CTj, 
where CTj's are Pauli matrices, on two-qubit systems, H' -<i^u H if and only if 
h' = [h'l, h'2, h = {hi, h2, h^) , where <s is the s-majorization defined in 

[11]. Thus we can see that in the case of efficient simulation of Hamiltonians on 
two-qubit systems, eigenvalues of Hamiltonians play a crucial role, since h and 
h! can be determined from the eigenvalues of Hamiltonians H and H' uniquely. 
However Theorem 1 implies that in the case of arbitrary bipartite quantum sys- 
tems, the algebraic-geometric invariants which are independent of eigenvalues 
play a more fundamental role. This is also illustrated in the following example 
of efficient simulation of Hamiltonians in H\ (g) . 

Example 1. Let ]wi), |t^2)|^^3) be the following 3 unit vectors in H\ H^. 



l^i) = 75(e'''Hll> + |22) + |33)) 
\v^) = ^^{e^^^\l2) + \2^) + (1) 
1^3) = 75(6^"^ |13) + |21) + |32)) 

,where 771, 772, % are 3 real parameters. Let Hfj^^j^^^rji = + \'U2){v2\ + 

|^^3)(^'3|- This is a continuous family of Hamiltonians in iJ^^^^^.^ya of rank 3 pa- 
rameterized by three real parameters. 

It is easy to calculate that (TJ^^ ,,,2 ,,,3 ) is just the elliptic curve (see [21] , [22] ) 

in defined by rf+r^+ rf - ';Z,'^f^llt,fJs rir2r3 = 0. Set 772, 773) = 

e.'lt^+rr+tr)"'? ^i^) = V^-H27)^^ ' ^^"^^ H9im,V2,V3)) is thc moduU function 
of elliptic curves. From algebraic-geometry it is known that if k{g(rii, 772, 773)) 7^ 
0,27, —216, then V^(-ffr)i ,772,773) is not the union of 3 lines and when k{g{rji,rj2, rjs)) = 
or 27 or —216, V^(-ff 771 ,772,773) is the union of 3 lines. Moreover V|(il,,i ,772,773) is 
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isomorphic to y|(ff^j^^^^^^) by projective linear transformations if and only if 
ki9{Vi,V2,V3)) = ^(S'WT'^L^s)) ( see section 7.2 , pp.363-396 of [22]). Thus we 
immediately know that Hq^.q cannot be efficiently simulated by -ffo.o,7r using 
local unitary operations from Theorem 1. Generally we have the following result. 

Corollary 1. iJ^j cannot be simulated by H^^^ ,^^ ,^^ efficiently by us- 
ing local unitary transformations,ie.,we cannot have Hyi^^^>^ ,^i^ <lu Hni,r)2,v3> 
kidivi^ 112,113)) f^idiv'i^V'iT'ns))! though the 3 nonzero eigenvalues of H^^^^^,^! 
Hr,[,n'2,ri'^ o-i^<i their partial traces are all 1. 

Proof. It is easy to calculate the eigenvalues to check the 2nd conclusion. 
The first conclusion is from Theorem 1 and the above-described well-known fact 
about elliptic curves. 

This example strongly suggests that the problem of efficient simulation of 
Hamiltonians on arbitrary bipartite quantum systems is quite different to the 
problem in two-qubit case as studied in [11]. 

Let S be the swap operator on the bipartite system H\ (g) ffg defined by 
= For any Hamiltonian H, S{H) = SHS^ corresponds to the Hamil- 
tonian evolution of H with A and B interchanged. It is very interesting to 
consider the problem if H can be simulated by S{H) efficiently . This led to 
some important consequences in the discussion VII of [11]. For example it was 
shown there are examples that H and S{H) cannot be simulated efficiently with 
one another in higher dimensions. Thus in higher dimensions nonlocal degrees 
of freedom of Hamiltonians cannot be characterized by quantities that are sym- 
metric with respect to A and B, such as eigenvalues. This conclusion is also 
obtained from the above Corollary 1. From Theorem 1 we have the following 
necessary condition about H -<lu S{H). 

Corollary 2. Let H be a semi-positive Hamiltonian on H\ ® Hg. Suppose 
H <LU S{H). Then V^{H) = V^{H) for k = 0, ...,n- 1. 

The following is a Hamiltonian ff on 3 x 3 system for which H cannot be 
simulated efficiently by S{H). 

Example 2. H = |<^i)(0i| + |</'2)(02| + where, 

|.^i) = ^(|ll) + |21) + |32)) 

l'^^^ = 7TW^I^'^+"l''^) (2) 
l'^3) = ;7=^(|13) + A|23)) 
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Then it is easy to compute that V^(-ff) is the sum of 3 lines in CP^ defined 
by n + r2 = 0,ri + vr2 = and n + Ar2 = for ^ A and both v, A axe not 1, 
and VK-ff) is the sum of 2 lines in CP^ defined by r2 = and rs = 0. Thus we 
cannot have H ~<Lir S{H). 

The following example shows that our results can lead to non-trivial con- 
straints without referring to eigenvalues even in the two-qubit case if local terms 
are not neglected (as in our definition). Although local terms may be physically 
irrelevant in the setting of [11] and [16], Example 3 below illustrates mathemat- 
ically how our results work. 

Example 3. Let H = + X2H2){H and H' = + 

-^2lV'2)(V'2l be two Hamiltonians on (g) ff^, where A's are any given positive 
real numbers such that trH = trH' and 

|Vi) = ^(|ll) + |22)) 
IV'2) = ;^(|11)-|22)) 
W) = ;^(|ll> + |22)) ^ ^ 
m = |12) 

Then we know H and H' are two rank 2 Hamiltonians. It is easy to compute 
that V^(i?) is the algebraic set of two points (1 : 0) and (0 : 1) in CP^ and 
V^{H') is the algebraic set of one point (0 : 1) in CP^. Hence we cannot have 
Lu H from the Theorem 1. 

We can now observe the compatibility of our necessary condition Theorem 
1 with the sufficient and necessary condition in two-qubit case in [11]. For two 
two-qubit Hamiltonians H and H' in their normal forms, i.e., H = TiihiCJi ® cji 
and H' — Y^ih[ai ai, it is proved in [11] that H' -<lu H if and only if 
h' = {h'i,h'2,h'^) h = (/ii, /i2, hs), i.e., h' is s-majorized by h. It is clear that 
the 4 eigenvectors of any two-qubit Hamiltonian in its normal form arc exactly 
4 Bell states. Thus if iJ = Tnhidi ® Oi is of the form T — ^^h, where T is 
a semi-positive two-qubit Hamiltonian, then the algebraic-geometric invariants 
of T are fixed, i.e., V^iT) is empty, and V^{T) is the set of 2 points when 
rank(T) = 2 and empty when rank{T) = 3 or 4. Thus we can see that our 
necessary condition Theorem 1 is void when applied to two-qubit Hamiltonians 
in their normal forms, the necessary condition in this paper is compatible with 
the main result in [11]. 

Actually the algebraic geometric invariants in [1] can be used to give more 
necessary conditions for the efficient simulation of Hamiltonians by using local 
unitary operations. 
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Theorem 2. Let H and H' he two semi-positive Hamitonians on H^®H^. 
Suppose that there exists a representation of H as a convex combination H = 

Yilqi\vi){vi\, with positive qi's and the Schmidt rank of \vi) is min{m,n}. More- 
over V^{H') is not empty. Then H' cannot be simulated by H efficiently by using 
local unitary operations, ie., we cannot have H' -<lu H. 

Proof. From the condition, there exist positive pi, ...,pt and local unitary 
operations Ui (g) Vi,..., Ut O 14, such that, ll\^^piUi ® ViH{Ui (g) Vi)^ = H' . It 
is clear that {U^ ® V{)H{U^ ® V,)'^ = ^^-^qJ\{U^ ® Vj)wj)((t/, ® V{)Vj\, and thus 

H' = S-fij=iftgi|(C/i ® Vi)vj){{Ui (g) Vi)vj\. From Lemma 1 in [20] as cited ( 
Lemma 1 ), range{H') is the linear span of vectors {Ui Vi)vj for i = 1, ...,t 
and j = 1, s. From the above description about the computation of V^{H'), 
we can compute it by choosing dim{range{H')) linear independent vectors in 
this set {{Ui (g) V-C)vx, [Ui (g) Vi)vs, {Ut (g Vt)vi, {Ut (g Vt)Vs}. Therefore 
we can choose one of these dim{range{H')) linear independent vectors to be 
{Ui ® Vi)vi, whose Schmidt rank is min{m,n}. From [1] and the definition , 
we know that V2{H') has to be the empty set. This is a contradiction and the 
conclusion is proved. 

Example 4. Let H = \v){v\ and H' = ^{\u-i){ui\ + \u2){u2\) be two Hamil- 
tonians on where 

H = -1^(111) + |22) + 133)) 
ki) = 7l(|ll) + |22)) (4) 
75(111) -|22)) 

It is clear that V^{H') is the set of one point (0 : : 1) in CP'\ thus 
nonempty. On the other hand H satisfies the condition in Theorem 2. Thus we 
cannot have H' -<lu H. 

For multipartite Hamiltonians on (g • • -(giJ™" , the definition H' -<lu H 
can be naturally extended as follows. We say that multipartite Hamiltonian H' 
can be simulated by H efficiently using local unitary operations, written as 
H' ~<Lu H, if there are positive real numbers pi, ■■■,Ps such that pi + .. -\-ps = 1 
and unitary operations U^ , ...,U} ,...,Ui , ...,U" on respectively, 
such that, H' = pi{Ul (g • • • (g Ul')H{Ul (g • • • (g U^)^ + ... +Ps{Ul (g • • • (g 
Ug)H{Ul (g • • • f/")^. Then we can use the algebraic-geometric invariants in 
[2] to give the following necessary condition. 

Theorem 3. Let H and H' be the semi-positive Hamiltonians on the mul- 
tipartite quantum system (g • • • (g iJ™" with the same rank. Suppose that 
Lu H, that is , H' can be simulated by H efficiently by using local unitary 
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operations. Then V^. ■ .4. (H) = V^. . .a. (H') for any possible k and any 
possible Ai^,...,Ai. , here the equality of algebraic sets means they are isomor- 
phic via projective linear transformations of the product of complex projective 
spaces. 



The proof is the same as the proof of Theorem 1. 



The following example is 3 qubit case. 



Example 5. Let H and H' be rank 4 Hamiltonians on (g) ig) Hq, 
H = |</.i)(<^i| + |<^2)(</>2| + |<^3)(<^3| + |.^4)(<^4| and H' = \cj>'M\ + \cl>'^){cj>'^\ + 
l'^3)('^3l + l'^4)('^4l, where, 



1^1) = ^(1010) -1011)) 
|'^2) = ^(|100)-|110)) 
l</'3) = 73(1001) -|101)) 
|</)4) = ^(1000) - |111)) 
l'^'i) = 7l(|000)-|100)) 
1-^2) = 73(1001) -|101)) 
l<A3> = 7l(|010)-|110)) 
|</)^) = ^(1011) - |111)) 

Then we can compute that Fj.^ (F) is the sum of CP^ x (1 : 0), (0 : 1) x CP^ 
and (1 : 0) X (0 : 1) in CP^ x CP\ and VX.g{H') is the set of two points 
(1 : 1) X (0 : 1) and (1 : 1) x (0 : 1) in CP^ x CP^. Thus from Theorem 3, we 
cannot have H' -<lu H. 



In conclusion, we have proved necessary conditions for the efficient simula- 
tion of both bipartite and multipartite Hamiltonians using local unitary oper- 
ations, which are independent of eigenvalues and based on algebraic-geometric 
invariants. These conditions indicated that the in higher dimension bipartite 
cases or multipartite cases, the relation of efficient simulation of Hamiltonians 
depends more on the eigenvectors than eigenvalues. This is quite different to 
the two-qubit case studied in [11]. 

It is natural to ask if the techniques from algebraic geometry used here can be 
extended to find not only necessary but also sufficient conditions of simulating 
semi-positive bipartite Hamiltonians by another with local unitary operations. 
For low rank semi-positive bipartite Hamiltonians or semi-positive Hamiltoni- 
ans on low dimensional bipartite systems, it seems that the eigenvalues of these 
Hamiltonians, eigenvalues of their partial traces and algebraic-geometric invari- 
ants are near a complete set of invariants under local unitary operations, ie., we 
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almost can determine exactly in which equivalent class the bipartite Hamilto- 
nians are if we know all these invariants. Thus in these cases it seems hopeful 
to extend the techniques here to find necessary and sufficient conditions of sim- 
ulation problem of Hamiltonians. However in general case we think it would 
be difficult to get necessary and sufficient conditions about this problem based 
on present-known invariants. We speculate that more invariants of bipartite 
Hamiltonians under local unitary operations have to be found for the purpose 
to completely describe the ability of bipartite Hamiltonians to simulate one an- 
other with the help of local unitary operations. 
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